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We investigate the fluctuations of the diquark-pair field and their effects on observables 
above the critical temperature T c in two-flavor color superconductivity (CSC) at moderate 
density using a Nambu-Jona-Lasinio-type effective model of QCD. Because of the strong- 
coupling nature of the dynamics, the fluctuations of the pair field develop a collective mode, 
which has a prominent strength even well above T c . We show that the collective mode is 
actually the soft mode of CSC. We examine the effects of the pair fluctuations on the specific 
heat and the quark spectrum for T above but close to T c . We find that the specific heat 
exhibits singular behavior because of the pair fluctuations, in accordance with the general 
, theory of second-order phase transitions. The quarks display a typical non-Fermi liquid 

behavior, owing to the coupling with the soft mode, leading to a pseudo-gap in the density 
' of states of the quarks in the vicinity of the critical point. Some experimental implications 

Ph! of the precursory phenomena are also discussed. 

§1. Introduction 

In recent years, there have been many studies of the properties of dense and cold 
matter. 1 ) In such a system, the baryonic density is so high that quarks and gluons 
are expected to be deconfined to make a quark matter. 2 ) Then, the attractive quark- 
quark interaction in some channels should give rise to a Cooper instability leading to 
color superconductivity (CSC) at sufficiently low temperatures. 3 ^' 4 ^ The asymptotic 
freedom of QCD allows us to use perturbation theory to describe CSC at extremely 
high densities with u > 10 8 MeV, where fi is the quark chemical potential. 5 -* In 
this region, the weak-coupling theory analogous to BCS theory, i.e. a mean- field 
approximation (MFA), is valid. In such a system composed of u, d and s quarks at 
extremely high density, it is believed that the quark matter takes the special form of 
the color-superconducting phase, i.e. the color-flavor locked (CFL) phase, 6 ) in which 
all kinds of quarks take part equally in the pairing. 
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What is the phenomenological significance of CSC including the CFL? Possible 
relevant systems consisting dense QCD matter include the cores of compact stars, 7 ) 
both in equilibrium and in the newly-born stage and the intermediated states created 
in the heavy-ion collisions that are expected to be performed in the forthcoming 
facilities at GSI and J-PARC. 8 )' 9 ) One of the basic points concerning real systems of 
this kind lies in the fact that the quark-number chemical potential /j, and hence the 
baryon density are moderate; fi should be at most ~500 MeV. Therefore the ideal 
situation realized in an extremely dense system may not be expected in these real 
systems but various complications come into play in the determination of the nature 
of CSC. 

To realize CSC in a compact star at vanishing temperature, the color- and 
electric-charge neutrality conditions must be satisfied as well as the /3-equilibrium 
condition. 10 )~ 12 ) These conditions, in turn, induce a mismatch of the Fermi mo- 
menta of quarks of different flavors and colors. This is the case in particular when 
the constituent mass of the strange quark, M s (which ranges from some 100 to 500 
MeV owing to the dynamical symmetry breaking of the chiral symmetry) is com- 
parable with ii. Thus we see that the incorporation of the charge neutrality and 
/^-equilibrium conditions, with finite M s taken into account, brings about quite in- 
teresting complications in the physics of CSC in quark matter at moderate densities. 
Indeed, it has been shown in the MFA that the combination of all these effects can 
lead to a variety of pairing patterns, including the so-called gapless CFL's. 13 )~ 20 ) 

In this paper, we focus on another important feature of QCD matter at mod- 
erate densities, i.e. the strong-coupling nature of QCD at low energy scales, which 
invalidates the MFA. This strong coupling may imply the significance of large fluctu- 
ations of the diquark-pair field, especially in the vicinity of the critical temperature, 
7^8), 21), 22) j£ - g no t ew0 rthy that some recent analyses of the RHIC experiments 
suggest that quasi-bound quark-antiquark states may be formed in moderately hot 
QCD matter. This reflects the strong-coupling nature of QCD, 23 ) although the pos- 
sible existence of hadronic modes above T c was suggested earlier in Refs. 24) and 25) 
*). 

It is thus natural to expect that heated quark matter at a moderate density 
may also accommodate pre-formed diquark pairs as a pre-critical phenomenon of 
CSC phase transition. Owing to the strong-coupling nature inducing large fluctu- 
ations of the order parameter, CSC in heated quark matter at moderate densities 
may have some of the same basic properties as the superconductivity in strongly 
correlated electron systems, such as high-T c superconductivity(HTSC), rather than 
usual superconductivity in metals. 8 )' 27 ) We notice that materials at T > T c show 
various types of non- Fermi liquid behavior, one of which is pseudogap formation, i.e. 
an anomalous depression of the density of states (DOS) N(lo) as a function of the 
fermion energy uj around the Fermi surface. 28 )~ 30 ) It would be intriguing to explore 
whether quarks in the heated quark matter near the critical point exhibit similar 
abnormal behavior. In this paper, we investigate pair fluctuations of CSC for T 
above T c and explore their effects on some physical quantities, including the quark 
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spectra in such quark matter near the critical point. 

The relevant quark matters we have in mind are those created by the heavy-ion 
collisions or those realized in the core of a proto-compact star just after a supernova 
explosion. In these systems, the temperature is considerably high, and for this reason 
the /^-equilibrium condition may not be completely satisfied, in contrast to the case 
in the interiors of cold compact stars. Therefore, the difference between the chemical 
potentials of the up- and down-quarks in these systems is smaller than that under the 
/3-equilibrium condition, and hence a two-flavor superconductor (2SC) will be favored 
over other pairings which incorporate strange quarks. Furthermore, the effect of the 
difference in chemical potentials induced by the neutrality conditions will be smaller 
at finite temperature than that at T = 0, because the Fermi surface is diffused for 
T ± O. 15 ) Thus, one may simply introduce the same chemical potential for all the 
kinds of quarks as a fair approximation to study CSC in such hot quark matter. 

There are two types of fluctuations in superconductors, i.e., the fermion-pair and 
gauge-field fluctuations, irrespective of whether they are electric or color supercon- 
ductors. Gauge-field fluctuations are known to make CSC phase transition first order 
in the weak coupling region. 4 )' 31 )~ 33 ) In this region, it is known that CSC is strong 
type-I superconductivity 31 ^' 34 ) in which the fluctuations of the gauge field dominate 
the pair-field fluctuations. On the other hand, as is argued in Ref. 34), CSC is ex- 
pected to be type-II superconductivity at lower density, i.e. for fi < 500MeV, *) 
where the fluctuations of the pair field dominate those of the gauge field in contrast 
to the weak-coupling case. In the present work, therefore, we simply ignore the gluon 
degrees of freedom, and examine the effects of pair fluctuations near T c , as done in 
Refs. 8) and 27): Because pair fluctuations are inherent in second- (or weak first-) 
order phase transitions, the results in the present work should hold irrespective of 
the different pairing patterns and the chemical potential combination, as long as the 
phase transition to CSC is second order or weak first order. We treat the system at 
relatively low density, where the strange quark degrees of freedom do not come into 
play; accordingly the pairing pattern is taken as 2SC throughout this paper. 

To investigate a system at relatively low temperature and density, a perturbative 
QCD calculation is inadequate, because of the strong coupling. Lattice Monte Carlo 
simulations for finite \i are still immature for the present purpose, although much 
progress has been being made in recent years. 35 ) Therefore, it is appropriate to adopt 
a low-energy effective theory of QCD. Because we are considering the situation in 
which the diquark pairing dominates over that of the gluons, an effective model 
composed solely of the quark fields may be adopted. Such effective models include 
the instanton induced model 36 ) and the Nambu-Jona-Lasinio (NJL) model. 37 )' 38 ) 
Note that the latter model is a simplified version of the former. Both models are in 
fact used to explore the phase structure at low density. 38 ) -42 ) One should note here 
that an NJL-type theory can be deduced as a low-energy effective theory for dense 
quark matter on the basis of the renormalization- group equations. 43 ) It is also known 

In Ref. 34), it is shown by extrapolation from the weak coupling region that CSC is type-II 
superconductivity when T c > 14 MeV for fx = 400 MeV. In our model, T c = 40 - 60 MeV for fj, = 400 
- 500 MeV, as shown in the next section, and thus one may conjecture that type-II superconductivity 
is realized. 
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that the phase structures obtained in these models are qualitatively consistent, and 
they are similar to those calculated using the Schwinger-Dyson equation with the 
one-gluon exchange interaction. 22 )' 44 ) In this work, we employ the NJL model to 
explore the fluctuations in CSC. 

In previous short communications, Refs. 8) and 27), we studied the precursory 
phenomena to CSC: We showed for a vanishing wave number that the dynamical 
fluctuations of the pair field have a prominent strength near the vanishing frequency 
up to T ~ 1.2T C , 8 ) and for this reason that these fluctuations can give rise to inter- 
esting precursory phenomena, such as a pseudogap in the quark DOS near T c . 27 )*) It 
is worth emphasizing that Ref. 27) is the first investigation to explore whether and 
how the quasi-particle properties of quarks are changed by the precursory paring 
mode in CSC. In the present paper, we present a formulation of the basic theory and 
give a detailed account of both analytic and numerical calculational procedures. We 
also present a more extensive study of the properties of the precursory pair field with 
finite wave-numbers and the mechanism through which the pseudogap appears. We 
investigate for the first time the effects of the precursory pair fluctuations on the spe- 
cific heat. It is found that the heated quark matter close to the critical point of CSC 
shows typical non-Fermi liquid behavior, owing to the precursory pair fluctuations, 
which form a soft mode. 

This paper is organized as follows. In §2, we introduce our model Lagrangian 
and present the phase diagram obtained in the MFA in the model. In §3, we inves- 
tigate the behavior of the pair fluctuations above T c using linear response theory. It 
is shown that the fluctuations of the pair field develop a collective mode with a large 
strength even well above T c . We show that the complex frequency of the collective 
pair-field moves toward the origin in the complex energy plane, which implies that 
the pair fluctuations form the soft mode of CSC phase transition. We calculate the 
spectral function of the pair fluctuations, p(k,u), as a function of the momentum 
k and energy uj in order to elucidate the spatial and temporal behavior of the pair 
fluctuations when the temperature is lowered toward T c . We also present the behav- 
ior of the dynamical structure factor. In §4, we discuss the effect of the soft mode 
on the specific heat c v above T c and show that c v increases when T is lowered to T c 
and eventually diverges at T = T c in accordance with the singular growth of pair 
fluctuations. We examine how the soft mode affects the single quark spectrum and 
changes the quasi-particle picture of the fermion above T c in §5. We start from a 
calculation of the single-quark Green function in the T-matrix approximation to in- 
corporate the effects of the pair fluctuations in the quark sector. It is shown that the 
pair fluctuations give rise to a large decay width for quarks near the Fermi surface, 
i.e. non-Fermi liquid behavior of the quarks near T c . It is further shown that the 
anomalous behavior leads to a pseudogap in the DOS of the quarks. The chemical 
potential dependences of the quark spectrum and the DOS are also examined. The 
final section is devoted to a summary and concluding remarks. 



*' Possible pseudogap formation in association with the chiral transition due to phase fluctua- 
tions was discussed previously. 45 ^ 
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§2. Model and Phase Diagram 

In this section, after introducing our model Lagrangian, we recapitulate the 
derivation of the thermodynamic potential and present the phase diagram obtained 
in the MFA. 41 ) This forms the basis for the succeeding investigation of the nature of 
the precursory pair fluctuations. 

We employ a Nambu-Jona-Lasinio (NJL) model with two flavors and three col- 
ors, as mentioned in the Introduction: 

C = i)i^ + C s + Cc- (2-1) 
Here the quark- antiquark and quark-quark interactions, Cs and Cc are given by 

C c = Gc{i>ilsT 2 \A^ C ){4> C ilsT 2 \Aip), (2-2) 

with (x) = C / ip T (x) and C = 27270 being the charge conjugation operator. The 
matrices t 2 and \a (A = 2, 5, 7) are the antisymmetric components of the Pauli and 
Gell-Mann matrices for the flavor SU(2)f and color SU(3) C , respectively. We take 
the chiral limit putting m u = = 0, since the properties of the diquark condensates 
are affected very little by the small quark masses. 

We choose the scalar coupling constant as Gs = 5.01 GeV™ 2 and the three- 
dimensional momentum cutoff A = 650MeV so as to reproduce the pion decay 
constant f w = 93 MeV and the chiral condensate {ipip) = (— 250MeV) 3 in the 
chiral limit. 46 -* There are several sources to determine the diquark coupling con- 
stant Gc, for example, the instanton- induced interaction and the diquark-quark 
picture of baryons. 47 -* The former gives Gc/Gs = 0.5, 42 ) while the values G c /G s = 
0.49-0.73 48 ) and G c /G s = 1.5 - 2 49 ) have been obtained using the latter model. In 
the present work we fix 0^/^5=0.62, i.e., G c = 3.11 GeV" 2 , following Refs. 8), 40) 
and 27), where Gc is chosen so as to reproduce phase diagram similar to that cal- 
culated in the instanton-induced interaction. 39 ) 

To determine the phase diagram, we have to first derive the thermodynamic 
potential fl = —T log Tr e~^ K with 

K = j (fx {^H7 • V)V> - M^rV -Cs- Cc] ■ (2-3) 

Here, Tr denotes a trace operation over the color, flavor and Dirac indices. To apply 
the MFA, we assume a finite diquark condensate for the 2SC pairing A and the 
quark-antiquark condensate M, 

A = -2G C (ip C r?p) , M = -2G S (^} , (2-4) 

with r = i^T 2 \2- Employing the MFA for Cc and Cs, the thermodynamic potential 
in the MFA per unit volume is given by 41 ) 



Qmf{M, A- T, m ) = ~ log Tre-P K 
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Fig. 1. The calculated phase diagram in the T-/i plane for our model. The solid and dashed curves 
denote the critical lines of a first- and second-order phase transitions, respectively. 



M 2 , |Z\| 2 A f d 3 p 



+ ^- 4 /(^^ + Tl0g ( 1 + e ^)( 1 + e ^ + ) 



+e_ + e+ + 2Tlog(l + e- /3e -) (l + e"^)}, (2-5) 

where 



Ep = Vp 2 + M 2 , t ± = E p ± l i 1 e± = ^l + \A\ 2 , (2-6) 

and f3 = 1/T. The thermodynamic potential £2mf realizes its absolute minimum as 
a function of M and A in the equilibrium state; accordingly the optimal values of 
M and A satisfy the stationary conditions 



dM 



df^MF 

L) and 



A dA 



= 0, (2-7) 

M 



which are actually self-consistent equations for the condensates and are called "the 
gap equations". 

Using Eqs. (2-5) and (2-7), we can determine the phase structure of the model. 
The phase diagram in the T-fi plane is shown in Fig. 1: The dashed (solid) curves 
denote the critical line for a second-(first-)order phase transition. The critical chemi- 
cal potential for the chiral-to-CSC transition at T = is ji = 316 MeV in our model. 
The second order normal-to-2SC phase transition occurs somewhere in the range 
T = 35 - 60 MeV for /x = 350 - 500 MeV. In the following sections, we explore the 
pair fluctuations and precursory phenomena above T c in this region of ji. 

For later convenience, we now explicitly write down the critical condition for the 
normal-to-CSC phase transition with vanishing chiral condensate (M = 0). We first 
give the explicit expression of the gap equation for A in Eq. (2-7): 



r d 3 P r i 



^ ■ -itanh^±), 
e+ 2 /' 



A = 8G C A / — ^ <^ — tanh -— + — tanh } , (2-8) 



where e± = y/\p ± | 2 + |Z\| 2 , because we have set M = 0. Dividing Eq. (2-8) by A 
and then setting A = 0, the condition for determining the critical temperature T c of 
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the normal-to-CSC phase transition is obtained as 

Q = ! - SGg [ ** I J- tanh ftfezjO + -i- tanh + ">} . (2-9) 
G 7 (2*)»\l>-/l 2 p + M 2 J 

with p c = 1/T C . Below we find that Eq. (2-9) plays a crucial role for understanding 
the anomalous behavior of the pair fluctuations near the critical point. 

§3. Pair fluctuations above T c 

In this section, we discuss the properties of the precursory fluctuations of the 
diquark-pair field in the normal phase on the basis of linear response theory. Because 
we limit our attention to the behavior of fluctuations in the normal phase, we set M = 
A = 0. It will be shown that a collective mode corresponding to pair fluctuations is 
developed near the critical point; we find that this mode is the soft mode for the CSC 
phase transition in the sense that the pole of the spectral function in the complex 
energy plane moves toward the origin as T is lowered toward T c . We then calculate 
the spectral function and the structural factor of the pair field with finite energy and 
momentum and examine their temperature dependence. We also discuss how the 
chemical potential [i affects the properties of the soft mode. 

3.1. Linear response of the pair field 

In this subsection, we apply linear response theory in order to investigate the 
diquark-pair fluctuations. 8 ** The presentation of the theory here is more formal than 
that given in Ref. 8). See also Ref. 8) for a more intuitive discussion of the linear 
response of the pair field. 

When we apply an external perturbation H ex (t) to a thermal equilibrium state 
at t = to, the expectation value of an arbitrary operator, (0(x, t)), deviates from the 
initial equilibrium value (0(x,t)) eq , where 0(x,t) = e lKt O(x)e~' lKt , with K being a 
time-independent operator. In the linear response theory, this deviation is given by 

(0(x,t)} - (0(x,t)) eq = i I ds([H ex (s),0(x,t)}) eq . (3-1) 

Here, ( ) eq represents the thermal expectation value without H ex (t). 
In order to study the pair fluctuations of CSC, we choose 

0{x) = -2G c ^{x)r^ c {x), (3-2) 

H ex (t) = f d 3 xA* ex (x, t) {^ c {x)r^{x)) + h.c, (3-3) 



where A* x (x,t) is a classical external source. Substituting these into Eq. (3T) and 

)eq 



using the fact that (tjjrip c ) eq vanishes in the normal phase, we obtain the diquark 



pair field induced by A* ex , 

A* nd (x,t) = -2G c (iP C (x,t)r^(x,t)), 

/OO f' 
dt' / d 3 x'D R (x,t;x',t')A* ex (x',t'). (3-4) 
-oo J 
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Here, the response function D R is given by 

D R (x, t; x', t') = -M(t - t'){[ijj(x, t)r^j c (x, t),^ c (x', t')riP(x', t')]) eq , 



I 



|^£«(fc,u;)e-- N'Je'M-"'), (3-5) 



where we have taken the limit to ~^ — oo. The Fourier transformation of Eq. (3-4) 
gives 

A* nd (k,oj) = -2G c D R (k,u)A* ex (k,u). (3-6) 

Equation (3-6) contains a significant amount of information regarding elementary 
excitations in the system. We first note that if the external field has the form 

A% c (x,t)=A* ex (k,u>)e- iut+ik ~ (3-7) 

then Eq. (3-6) implies that there appears an induced pair field 

A* nd (x,t) = A* nd (k,co)e- M -* (3-8) 

with amplitude A* nd (k,uj). If the system has an intrinsic collective excitation with 
the dispersion relation w = uj s (k), the induced pair field will have a larger ampli- 
tude A* nd (k,uj) for uj ~ uj s (k) owing to a resonance mechanism; this implies that 
D R (k, uj) increases as uj approaches uj s (k) on account of Eq. (3-6) and eventually may 
diverge at uj = uj s (k). Conversely, if the amplitude A* nd (k, to) has a finite value with 
an infinitesimally-small external disturbance A* x {k,ui), the system would have an 
intrinsic collective excitation with frequency uj and wave number k. This situation 
is realized when the response function D R (k,co) is divergent, which may occur for a 
special uj with a given k. Thus we find that the equation 

D R (k,u J y 1 = (3-9) 

gives the dispersion relation uj = uj s (k) of the possible collective excitation in the 
diquark channel of the system. A remark is in order. Note that the frequency uj s (k) 
with given k can be a complex number; for instance, (i) if Re uj s (k) ^ 0, the mode is 
oscillatory, and (ii) if u> s (k) is pure-imaginary, the mode is diffusive. Below we find 
that Re oj s (k) is finite, but so small that the mode is almost diffusive. 

The information concerning the strength of the pair fluctuations is contained 
in the spectral function p(k,uj), which is defined through the response function 
D R (k,uj): 

p(k,uj) = --lmD R (k,uj). (3-10) 

7T 

The strength of the fluctuations is also represented by the dynamical structure 
factor, defined by 

= ~ ^TTi^) Im DR ( k ^), (3-11) 
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Fig. 2. The diquark propagator in the random phase approximation. 

which is positive definite. At finite temperature, S{k,uj) can be directly observed in 
scattering experiments; 50 ^ the van Hove scattering formula 50 ^ gives the cross section 
in terms of the dynamical structure factor. As we see below, p(k,u) represents the 
total width of the fluctuations, i.e., the difference between the decay rate and the 
production rate, and it is negative in the case to < for bosonic excitations. Accord- 
ingly, S(k,cj), not p{k,uj), represents the excitation probability of the fluctuations 
of the system at finite temperature. 

To calculate D R (k,uj), we employ the imaginary time formalism; we first calcu- 
late the response function in the imaginary time formalism (the two-particle Mat- 
subara Green function), 

2?(aj,r;0,0) = -{T T $(x, r)r^ c (x, r)^ c (Q, 0)/>(0, 0)) 

JLSLvfauJe^e** (3-12) 

with v n = 2-Kn/P being the Matsubara frequency for bosons. 

In order to evaluate Eq. (3-12), we employ the random phase approximation 
(RPA) in which we sum up the ring diagrams shown in Fig. 2, with the lines repre- 
senting the free Matsubara Green function in the normal phase (M = A = 0), 

Go(k,u; n ) = 1 — - , (3-13) 

(lLU n + /i)7o - k ■ 7 

and the vertices correspond to Cc- In this approximation, the Matsubara Green 
function (3-12) is given by 

V(k, u n ) = - y ' — -, 3-14 
2 1 + G c Q(k,v n ) 

where Q(k, v n ) is the contribution from the one-loop particle-particle correlation 
function, 

Q(k, v n ) = -2tJ2 I ?S Tr [crg (k - P ,u n - Lu m )rcgT( P ,Lu m )} . 

(3-15) 

Here ui n = (2n + l)vr//3 denotes the Matsubara frequency for fermions. The response 
function in the real time Eq. (3-5) is given by the analytic continuation iv n — > u + irj, 
and we obtain 

D R (k,u)=V(k,u n )\. M . = \ - r , (3-16) 
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with Q R (k,oj) = Q(k, v n )\iv„=u>+iri- It is thus seen that the equation to determine 
the dispersion relation for the collective excitation given in Eq. (3-9) is reduced to 

l + G c Q R (k,oj) = 0. (3-17) 

With some manipulations (see Appendix A for details), we find the following 
simplified form of Q R {k,oo): 



Q R (k,u) = N f (N c -l) J 



d 3 p 1 

(27r) 3 eie 2 



x i [(«, + „,» - ^ i !-/-(■»)-/-(«.) + ' - t+M - r+M 



oj + 2fi - ei - e 2 + if] oj + 2/x + e\ + e 2 + ir? 

_ [(ei _ _ A ( n«)-rto_ + f+M-rM \ ) . 

\cj + 2/i + e\ - e 2 + if] lu + 2fi - ei + e 2 + if] J ) 

(3-18) 

Here / ± (x) = 1/ [e^ 1 * 1 ^) + 1] is the Fermi-Dirac distribution function for the quark 
and anti-quark, respectively, and e\ = \p\, e 2 = p\ and k = |fe|. *^ The imaginary 
part of Q R (k,co) is given by 

x [(1 - /- ( ei ) - /-(e 2 )) <y( w + 2 M - ei - e 2 ) 

- (1 - /+( ei ) - /+(e 2 )) <5(o; + 2^ + ^ + e 2 ) 
+ (.T (ei) - / + (e 2 )) <5(a; + 2 M - ei + e 2 ) 

- (/ + (ei) - /" (e 2 )) S(lo + 2 M + ei - e 2 )] . (3-19) 

Each term in the square brackets of Eq. (3T9) corresponds to the decay process 
shown in Figs. 3 (a)-(d) and its inverse. To see this, the following identities for 
arbitrary quantities (i = 1,2) are useful: 1 - /i - / 2 = (1 - /i)(l - / 2 ) - /i/ 2 
and fi — / 2 = (1 — / 2 )/i — / 2 (1 — fi). In general, the imaginary part of the retarded 
Green function is related to the net decay rate or the total width. Because of the 
delta function in each term owing to energy-momentum conservation, the ranges of 
energy and momentum over which each decay process can take place are restricted to 
the region shown in the right panel of Fig. 3; for example, the decay process into two 
particles, as shown in (a) of Fig. 3, and its inverse occur only for oj > \k\ — 2/j,. The 
processes (b) and (c), the so-called Landau damping processes, and their inverses 
take place only for —\k\ — 2fi < oj < \k\ — 2[i\ hence the collective modes with 
vanishing momenta do not decay through these processes. We now see that the 
decay process (a) plays the dominant role near T c , because the strength of the pair 
fluctuations is concentrated near (uj, k) = (0, 0) for T ~ T c , as we show in the next 
subsection. We note that D R (k,u>) is analytic around (oj, k) = (0,0) for T > T c , 
because the functions Q R (k,oj), and hence D R (k,oj), have discontinuities only on 
the boundaries oj = ±\k\ - 2fi. **) 

Using the change of variables p <-> k — p, Eq. (3-18) is converted into Q R (k,oj) of Rcf. 8). 
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(a) 



(b) 



(c) 



(d) 




Fig. 3. The left panel shows the kinetic processes contained in \mQ R {k,ui)\ the corresponding 
inverse processes are not shown. The dashed line denotes the diquark pair field and the solid 
line the quark or the antiquark. The right panel represents the energy-momentum regions in 
which the decay processes (a)-(d) (and their inverse ones) occur. 



Here, we would like to discuss the critical behavior of the response function 
Eq. (3-16). In general, when the temperature approaches T c of the second-order 
transition from above, the fluctuations of the order parameter with a low frequency 
(small uj) and a long wave-length (small k) become easily excited. This implies, on 
account of Eq. (3-6), that D R (k, uj) with small uj and k becomes larger as the system 
approaches the critical point. At T c , the system becomes unstable with respect to 
uniform diquark-pair formation, and a finite and permanent pair-field is formed with 
an infinitesimally-small external field. Then, according to Eq.(3-6), D R (0,0) should 
be divergent at T = T c . One can show that this is indeed the case, that is 



D R (0,0) 



-i 



T=T C 



= 0. 



(3-20) 



In fact, the denominator of Eq. (3-16) for vanishing k and uj becomes 
1 + G C Q R (0,0) 



8G C 



(27T) 



^tanhfc^ + ^tanh^) 



p- \i 



p + H 



(3-21) 



which is found to vanish at T = T c , on account of Eq. (2-9). Equation (3-20) is called 
the Thouless criterion, 51 ) which may be used to determine the critical point. 

For the momentum integration of Eq. (3-18), we employ the following cutoff 
scheme. 52 )> 53 ) First, we note that the imaginary part of Q R (k,uj) is free from an 
ultraviolet divergence. Thus, we can evaluate the imaginary part without introducing 
a cutoff. In this scheme, the imaginary part, Eq. (3-19), can be reduced to a compact 
form (see Appendix A), 



Im Q R {k,uj) = —Nf(N c - 1)T 



{uj + 2fi) 2 - k 2 
2irk 



Equation (3-18) is more convenient for the following calculation than that given in Ref. 8). 

**' In the case of particle-hole modes, a discontinuity of the response function exists at u> 
±|fc|. 37 ' Hence, the response function is not analytic at the origin. 
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coshfoj + k)/4T „. . , w 1 

Then, we evaluate the real part of Q R (k,uj) from the imaginary part by using the 
dispersion relation. We introduce a cutoff at this stage, because Re Q R (k,u>) has 
an ultraviolet divergence. The cutoff should be chosen so as to satisfy the Thouless 
criterion for k = 0, discussed above. Thus, as shown in Appendix A, the real part 
of Q R is expressed as 

R e Q«( fc , W ) = -ip ^ QR{k fK (3.23) 

where P denotes the principal value.*) This cutoff scheme has the advantages that 
the imaginary part of Q R (k,u>) does not destroy any conservation laws through the 
introduction of a cutoff and it does reflects the symmetries of the system. In addition, 
the simple form Eq. (3-22), which is derived without introducing a cutoff, is quite 
convenient for numerical calculations. 

3.2. Collective excitation 

In this subsection, we obtain the dispersion relation of the collective mode in the 
diquark channel by solving Eq. (3T7) for to, with k given. Solutions of Eq. (3T7) 
should exist in the lower-half complex energy plane C~, because the imaginary part 
of the pole should be negative; otherwise the system would be unstable with respect 
of the creation of collective modes. Thus we actually solve the equation, 

1 + G c Q R (k,z) = (3-24) 

for a complex variable uj = z with Q R (k,z) = Q R (k,co)\ ul ^ z . In fact, it can be 
numerically verified that Eq. (3-24) does not have any solution in the upper-half 
plane C + . In order to calculate Eq. (3-24) for z G C~, we must perform the analytic 
continuation of Q R (k,Lo) to C~; a simple analytic continuation of Eq. (3T5), defined 
as 



Q(k,z) = Q{k,u n )\ Wn ^ z , zGC, (3-25) 

with C denoting the entire complex plane, has a cut along the real axis. Hence 
Q R (k,z) ^ Q(k, z) in C~. The retarded function Q R (k,z) in C~ is defined on an- 
other Riemann sheet of Q(k, z), and therefore one has to derive the analytic function 
of Q(k,z) in this Riemann sheet to solve Eq. (3-24). As shown in Appendix B, the 
form of Q R (k, z) in C~ is found to be 

Q R (k,z) = Q(k,z) + 2il(k,z), (3-26) 

*^ The Thouless criterion does not restrict the range of integration in Eq. (3-23) for k 0, 
because this criterion is a property of static homogeneous matter. In this sense, we have no criterion 
to determine the range of integration in Eq. (3-23) for k 7^ 0. In this work, however, we simply 
assume that the range does not change even for finite momentum. In any case, the final results for 
the relevant range of |fc| and ui do not depend on the choice of cutoff scheme. 
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Fig. 4. The pole u) s of the collective mode above T c for vanishing momentum with various p. The 
points denote the positions of the poles at the reduced temperature e = (T — T c )/T c = 0, 0.1, 
0.2 ■ ■ • . Note here that u> s (k = 0) | e =o = 0. The pole approaches the origin as e is lowered. 



with I(k,z) = lmQ R (k,u)\ u -> z . 

Using Eq. (3-26), we find a solution of Eq. (3-24) in C~. Note that the Thouless 
criterion Eq.(3-21) ensures that there exists a solution at the origin for T = T c , 
irrespective of \i. The temperature dependence of the solutions with k = is already 
presented in Ref. 8). We display them in Fig. 4 for three chemical potentials, \i = 350, 
400 and 500 MeV; the three curves correspond to these three values of fx, while the 
dots on the curves denote the position of the poles for each reduced temperature, 
e = (T — T c )/T c = 0, 0.1, 0.2, • • • . It can be seen that the pole goes up to the origin 
in the complex energy plane as T is lowered toward T c . It is also seen that the rate 
at which the origin is approached is almost the same for all the chemical potentials. 
The result clearly shows that the collective mode of the pair fluctuations is the soft 
mode 37 ^ of CSC phase transition, which gives rise to the precursory phenomena 
of CSC phase transition, as shown in the following sections. We will see that the 
softening of the pole causes a rapid growth of the peak in p(k,uj) and S(k,u) as 
T^T C . 

It is noteworthy that Re uj is small but finite. In the weak coupling limit, as 
described by BCS theory, the pole of the soft mode appears on the imaginary axis, 
i.e., Re uj = 0. This is due to a particle-hole symmetry inherent in the weak coupling 
limit: In metal superconductors, an attractive interaction exists only for fermions 
with energies satisfying \uj\ < ujd, as measured from the Fermi energy Ep with the 
Debye frequency ojd, which satisfies up, <C Ep. Thus, the DOS for the fermions 
participating in the pairing can be treated as a constant in this energy range, and 
accordingly the range is symmetric with respect to the Fermi surface. For CSC, the 
pairing involves all of the states in the Fermi sphere, not just those around the Fermi 
surface. This leads to an asymmetry between particles and holes. This asymmetry 
is the origin of the finite real part of the soft mode. From Fig. 4, we find that the 
ratio Re oo/Im uj becomes smaller for larger fi. This is plausible, because the rate 
of the particle-hole asymmetry becomes smaller as \i is increased. It is instructive 
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k [MeV] Re co s [MeV] 



Fig. 5. Momentum dependence of the pole co B (k) for /i — 400 MeV and s = 0, 0.2 and 0.4. The left 
panel depicts the k dependence of Imu> s (fc) . The right panel depicts the momentum dependence 
of uj s (k) in the complex energy plane. The pole moves away from the origin as k is increased. 



to note that an increase of Im u, along with that of n, is also understood in a 
different way: The density of states on the Fermi surface into which the collective 
mode decays increases with the chemical potential. Thus, a larger decay probability 
of the collective mode is realized, and this leads to a larger Im uj. 

The fact that the pole has both real and imaginary parts implies that the dy- 
namical behavior of the order parameter near T c is a damped-oscillator mode, while 
the pure imaginary poles in the weak coupling limit correspond to over-damped 
modes, or diffusive modes as mentioned above. However, the absolute value of Im u 
is larger than that of Re co even in the present case. This implies that the dynamical 
behavior of the order parameter of CSC is approximately the same as that of the 
over-damped case, and it is accurately described by a non-linear diffusion equation 
like the time-dependent Ginzburg-Landau (TDGL) equation in the weak coupling 
theory. 8 )' 54 ) 

Next, let us consider the momentum dependence of the pole oj s (k). In Fig. 5, 
we plot Imu} s (k) as a function of k. It is seen that |Imu; s (A;)| becomes larger as k is 
increased. This means that the lifetime of the fluctuations of the pair field becomes 
shorter with smaller wavelength. The pole u: s {k) in the complex energy plane is 
shown in the right panel of Fig. 5. It can be seen that the ratio Rew/Imu; is almost 
independent of k and e. 

3.3. Softening of the pair fluctuations 

In the previous subsection, we saw that quark matter at values of T above 
but near T c develops a collective mode owing to diquark-pair fluctuations. In this 
subsection, we show how the strength of the fluctuations changes when T is lowered 
toward T c . 

We first show the temperature dependence of p(k, to) at vanishing momentum 
k = for three different chemical potentials, fi = 350, 400, 500 MeV, and at several 
reduced temperatures in Fig. 6. 8 ) One can see that the position of the peak moves 
toward the origin, while the peak height increases as the temperature is lowered 
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Fig. 6. The spectral function p(k,cu) with k = for /i = 350, 400 and 500 MeV and several values 
of e. The height of the peak increases as T is lowered in the normal phase for all values of fi. A 
clear peak is seen even at as high a temperature as e — 0.2. 



toward T c , although the rate of growth of the peak decreases as /i is increased. 

Using the Thouless criterion Eq. (3-21), it can be analytically shown that the 
spectral function has a divergent peak at the origin, (k,co) = (0,0), for T = T c . It 
can be seen in Fig. 6 that the peak remains quite distinct as away from T c as e = 0.2. 
In the case of electric superconductivity in metals, the effect of fluctuations is small, 
because of weak coupling. Precursory phenomena in electric conductivity can be 
observed experimentally. For example, there appears anomalous enhancement called 
'paraconductivity' in the electric conductivity just above the critical temperature. 
However, the range of T over which clear paraconductivity is seen is limited to those 
satisfying e < 10~ 3 , even in superconductors with large fluctuations such as dirty 
alloys and low-dimensional materials. Therefore, we can conclude that fluctuations 
in CSC survive for values of e in two or three orders larger than in the case of electric 
superconductors. It follows that the precursory fluctuation phenomena of CSC can 
exist over a rather wide range of T, and hence they can be used as experimental 
signatures of CSC. 

The above arguments are based on the behavior of the spectral function p(k, to). 
We can draw further conclusions on the strength of the fluctuations if we consider 
the dynamical structure factor S(k,ui), defined in Eq. (3T1). As mentioned above, 
S(k,uj) is more directly related to observables than the spectral function p(k,co), 
which has a negative value for bosons when u < 0, as shown in Fig. 7. We display 
the s dependence of S(k, u) for p = 400 MeV at vanishing momentum transfer in 
Fig. 7 (b). It can be seen that a peak appears in S(k,u>) and moves toward the 
origin as the temperature approaches T c . This behavior is consistent with that of 
the spectral function, and it is a reflection of the softening of the pair fluctuations. 
We remark that the position of the pole corresponds to the peak of the dynamical 
structure factor, not the spectral function. 

Next, we show the energy-momentum dependence of S(k, u) at (i = 400 MeV 
and e = 0.02 and 0.1 in Fig. 8. The peak near u = flattens as k is increased in 
each figure, which means that the fluctuations with smaller wavelength tend to be 
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Fig. 8. The energy and momentum dependence of S(k, u>) for fi — 400 MeV. Here, we have e = 0.02 
(e = 0.1) in the left (right) panel. The peak in S(k,u>) grows, which implies that the pair 
fluctuations around to = k = become stronger, as e approaches 0. 



suppressed. 

§4. Specific heat 

In the previous section, we showed that there exist strong pair fluctuations even 
well above T c . In the present and following sections, we evaluate the effects of the 
pair fluctuations on some observables. In this section, we calculate the specific heat, 
taking into account the precursory fluctuations above T c , and show that the pair 
fluctuations cause an anomalous enhancement of the specific heat over range of T 
that is similar to that over which a prominent peak in the spectral function is seen. 

The specific heat per unit volume, c v , is calculated from the thermodynamic 
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^0= Q, 

An = {yy +(^yy + (v^i + ••• 

Fig. 9. The diagrams taken into account in the calculation of the thermodynamic potential J? = 
J?o + Ofi,. Here, J?o is the thermodynamic potential of a free fermionic system, and Qr, is the 
contribution from the fluctuations of the pair field. 



potential per unit volume, Q, as 



= ~T^. (4-1) 



Therefore, we first evaluate the thermodynamic potential I? in the normal phase, 
incorporating the pair fluctuations: 

f2 = f2 + {2 a .. (4-2) 

Here, Qq (i?fl.) denotes the contribution at the mean-field level (from the pair fluctu- 
ations). Then, the specific heat c v is also divided into the two parts, corresponding 
to i?o and f2&, respectively: 

Cv = C ° + C fl-. (4-3) 
The thermodynamic potential of the free quarks Qq reads 



[2 = Q M f{M = 0,A = 0) 
rfik 

—-IHog&CMn). (4-4) 



Here, f2± consists of the summation of the ring diagrams shown in Fig. 9, where all 
vertices correspond to the diquark interaction term Eq. (2-2). Below, we find that 
these diagrams correspond to those considered in the response function D R (k,uj) 
in the RPA, and hence fi^ exhibits anomalous behavior through D R (k,uj). The 
lowest-order diagram in is calculated to be 

Q { i ] = -2G C T 2 Y, I E Ti[i l5 T 2 X A g(p 1 ,Lo m )i 1 ,T 2 X A g(p 2 ,u; n )] 

m,n J ' A=2,5,7 

= -2G c T 2 Y, I ^^tj 2 3Tr ^5r 2 A 2 g(p 1 ,a; m )H 5 r 2 A 2 g(p 2 ,a; n )] 

m,n ^ ' 

= 3T Ej^ G cQ(P^n). (4-5) 



In the second equality here, we have used the fact that the A = 2, 5 and 7 terms in 
the first line of Eq. (4-5) give the same contribution; we have incorporated them all 
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Fig. 10. The contour of integration in Eq. (4-8). 



into the A = 2 term with an overall factor of 3. Physically, this factor corresponds to 
the existence of three degenerate collective excitations of the pair field in the normal 
phase, that is, the red-blue, blue-green and green-red collective modes in the color 
space. This factor appears in the diagrams of J7fj. at all orders. Similarly, the Z-th 
order diagrams of i? n . can be evaluated as 

^ = -J T Ej^i- G cQ(P^)] 1 . (4-6) 
Summing up all the terms, we obtain 

= E * = " 3T E / T^E! 1-GcQ<P, »n)] 1 (4-7) 

1=1 n J > 1=1 

= 3T E/ ^Iog[l + Gcfi(p,^)] 

= ! / (^ P / ^^TT Iml °g + ^(P.o,)] • (4-8) 

In the last equality, the contour of integration is modified so as to avoid the cut along 
the real axis, as shown in Fig. 10. 55 ) Note that the argument of the logarithmic 
function is the same as the denominator of the response function of the pair field 
D R (k,uj), which, we have seen, exhibits a singular behavior near T c . Therefore, i?fl. 
should also show an anomalous behavior near T c , which gives rise to an enhancement 
of C A-. 51 ) 

Before taking the derivatives in Eq. (4-1), we expand the argument of the log- 
arithmic function in Eq. (4-8) in a Taylor expansion about u> = k = and T = T c . 
To lowest order, this yields 



Gc 1 + Q R (k,oj) ~ C cj + C\k\ 2 + Ae, 



(4-9) 
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Fig. 11. The specific heat per unit volume in the free fermionic system, c°, and the effects of the 
fluctuations of the pair field, c^'. The total specific heat of the system is c v = c° + c^'. It is 
seen that as e decreases, the enhancement of c^' begins to become appreciable near e — 0.05-0.1. 
The dot-dashed curve represents the specific heat obtained with only the static part of fin.- 



with 

,dQ R (0,0) 



n OQ R (0,0) 

On = 



c ^ 9Q^(0,0) 



T=T C 



d\k\ 



, A=T- 



/ 01' 



. (4-10) 



T=T C 



Note that there is no zeroth-order term in the expansion Eq. (4-9). This follows 
from the Thouless criterion, Gq 1 + Q R (0, 0)| t=t c = 0. One can numerically check 
that the right-hand side of Eq. (4-9) indeed does provide a good approximation of 
the behavior of the left hand side over a rather wide region of to, k and e; roughly 
for those values satisfying \co\, \k\ < 120 MeV and e < O.3. 56 ) As was shown in the 
previous section, the strength of the pair field is concentrated around to = \k\ = 
as T approaches T c . Therefore, the fluctuating modes in the vicinity of the origin of 
the uj-k plane are expected to give the dominant contribution to the thermodynamic 
potential and hence the specific heat. Thus, use of the simple expansion given in 
Eq. (4-9) in their calculation is justified. 

In Fig. 11, we plot the behavior of the specific heat c° and c^', above T c . It is 
seen that c^' diverges as T — > T c . A clear enhancement of c^' is seen already near 
£ = 0.05-0.1. The range of temperatures in which is larger than c° is called the 
Ginzburg-Levanyuk region. 57 )> 58 ) In our case, the Ginzburg-Levanyuk region exists 
up to e = 0.02 above T c . In electric superconductors, the Ginzburg-Levanyuk region 
is continued only to values e satisfying e < 10~ 3 , even for dirty alloys. Therefore, 
the Ginzburg-Levanyuk region for CSC is very much larger than that for electric 
superconductors. It can be numerically checked that c^' ~ e -1 / 2 near T c . This 
critical exponent of c v is the same as that obtained from the Ginzburg-Landau 
equation without nonlinear terms. 59 ) This is because we have adopted the RPA 
for the pair field. However, the RPA is not valid in the vicinity of T c , because there 



20 



M. Kitazawa, T. Koide, T. Kunihiro and Y. Nemoto 



the amplitude of the fluctuations becomes large, and therefore the nonlinear effects 
begin to play a significant role. For this reason, the true critical exponent of the 
specific heat should differ from —1/2. 

The singular behavior of the specific heat above T c of the superconductivity is 
usually studied using only the static part of .Ofl., 59 ) i.e., the n = component of the 
Matsubara sum in Eq. (4-7). However, we have explicitly included the dynamical 
modes as well. We now give an account of the relation between the two approaches. 
The most singular behavior of i?fl. which causes the enhancement of c^' near T c 
certainly comes from the static part of i?fj., because only the static part in Eq. (4-7) 
diverges as e — > 0. Therefore, the contribution from the static part dominates the 
behavior of Or., at least near T c . In Fig. 11, the specific heat obtained using only the 
static part is plotted by the dot-dashed curve. It is seen that the singular behavior 
of near T c is reproduced by the contribution from the static part of alone, 
and the contribution from the n^O components of i?fl. to is smaller than c°. 

Recently, Voskresensky calculated the specific heat above the CFL phase and 
found that the Ginzburg-Levanyuk region exists up to e ~ l, 9 ' which is more than 
one order of magnitude larger in the units of e than that found in the present case. 
The origin of this difference can be understood as follows. In two-flavor quark matter, 
the degeneracy of the collective modes gives rise to the factor of 3 in Eq. (4-8), while 
the existence of nine collective modes above T c for the CFL phase can cause a factor 
of 9 to appear in c^, provided that the strange quark mass is ignored. Therefore, 
above the CFL phase can take a value three times larger than that in ours. The 
coefficients of the Ginzburg-Landau equation in Ref. 9), which are determined in 
the weak coupling limit, further increase c^'. These factors can account for the wide 
Ginzburg-Levanyuk region obtained in Ref. 9). 

§5. Non-Fermi liquid behavior due to fluctuations 

As seen in the previous sections, the pair fluctuations form a well-developed 
collective mode near T c in CSC, and they cause anomalous behavior in the specific 
heat. In this section, we explore how the pair fluctuations in turn affect the properties 
of quarks. We show that the pair fluctuations give rise to non-Fermi liquid behavior 
of the quarks near T c : It is shown, for instance, that a large decay width is acquired 
by the quarks near the Fermi surface. It is also shown that the anomalous behavior 
leads to the pseudogap in the DOS of the quarks. 

The pair field may exhibit both amplitude and phase fluctuations. We assume 
that the amplitude fluctuations dominate the phase fluctuations. 29 ) Thus we are led 
to employ the T-matrix approximation, 60 ) which is suitable to evaluate the effects 
of the amplitude fluctuations of the pair field on the quarks. Historically, the DOS 
above T c was first calculated within the T-matrix approximation for the case of 
the electric superconductivity in the weak coupling regime more than three decades 
ago. 61 ) -63 -* It was shown that a pseudogap can form but only in the vicinity of T c . 
After the discovery of the HTSC, the same approximation was reconsidered in the 
strong coupling regime in the study of the DOS for HTSC. 30 )' 64 ) In this section, 
we compute the quark Green function in the T-matrix approximation and evaluate 
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the spectral function and the dispersion relation of the quarks in the relativistic 
kinematics. 

5.1. T-matrix approach 

The one-particle quark spectral function is defined by 

A(k,u) = -- -lmG R (k,uj) 

7T 

1 G R (k,Lo)-G A {k,io) 1 G R (k,m) -~/°G R Hk,uj)"f 

~ 7T 2Z " 7T 2i ' (54) 

where G R (k,uj) and G A (k,u>) are the retarded and advanced Green functions of the 
quarks, respectively. Conversely, G R (k,oo) is given in terms of A(k,u>) as 

G R (k,oj) = [ dJ Stl^l , ( 5 . 2 ) 
J oj — oj + irj 

The spectral function A(k, to) has a Dirac matrix structure in the relativistic formal- 
ism. From the rotational and parity invariances, the Dirac indices of A(k, uo) can be 
decomposed into three parts, 

A(k, uj) = p (k, w)7° - p v (k, Lo)k ■ 7 + p s (k, oj), (5-3) 

with k = k/\k\. Actually, p s (k,oj) vanishes in the chiral limit, which we have taken. 
Note that po(k,oj) is related to the number density of quarks, and hence to the DOS 
of the quarks, as 

/d s k 
—Tr c ,f[p (k,u)}, (5-4) 

with Tr c f denoting the trace over color and flavor indices. 

To calculate the quark Green function Eq. (5-2), we use the Matsubara for- 
malism. The Dyson-Schwinger equation for the quark Matsubara Green function 
G(k, oj n ) reads 

Q(k,oo n ) = g (k,oj n ) + g (k,oj n )IJ(k,oj n )g(k,oj n ), (5-5) 

where U(k,oj n ) is the quark self-energy in the imaginary-time formalism. We take 
the following approximate form for the self-energy ZJ(k,oj n ) 



— 3 S{p + k, oj n + «/ m )Go(k, oj'J (5-6) 



(2ir) 

where the T-matrix S(k, u n ) is 



o9k ~ 

■■—~(p + k,oj n + oj'jg (k,oj' m ), (5-7) 



S(k, v n ) = - G C l + GcQ{k Vn) = ~Gc (1 - 2G c V{k, u n )) , (5-8) 
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Fig. 12. The Feynman diagrams representing the quark Green function in the non-self-consistent 
T-matrix approximation. The thin lines represent the free propagator Go, while the bold lines 
represent the full propagator Q. 

with T>(k, u n ) defined in Eq. (3-14). This is the so-called non-self-consistent T-matrix 
approximation, which is consistent with RPA employed in §3 to describe the pair 
fluctuations. Figure 12 is the schematic representation of our T-matrix approxi- 
mation; the thin and bold lines represent the free and full Matsubara propagators, 
Qo(k,u>) and Q(k,oo). If the thin lines in the diagram are replaced by the bold 
lines, the approximation becomes the self- consistent approximation. We comment 
on this approximation in the concluding remarks. The T-matrix S(k, v n ) show the 
same anomalous behavior as the response function near T c , in accordance with the 
softening of the pair fluctuations. 

In the second equality of Eq. (5-7), we have used the fact that the sum of the 
Gell-Mann matrices in Eq. (5-6) yields ^2a(Xa) 2 = 2. This factor corresponds to the 
two possible patterns of E(k,uj n ) in the color space owing to the existence of three 
degenerate collective modes of the pair field in the normal phase; a red quark, for 
example, can create a red- green and red-blue collective modes incorporating green or 
blue quarks, respectively. In the case of the specific heat discussed in §4, the factor 
3(= N c ) appears because the three degenerate collective modes contribute equally. 
On the other hand, only two of the three collective modes can couple to an incoming 
quark, and this is the reason that the factor 2(= N c — 1) arises in the self-energy 
appearing in Eq. (5-7). 

The analytic continuation of S(k, oj n ) to the real axis from the upper-half com- 
plex energy plane gives the self-energy in real time, 



After several manipulations (summarized in Appendix C), this self-energy is found 
to be 



E R (k,u;) = E(k,u n )\ ii 



(5-9) 




(5-10) 



with 



= E(k,v n )\ 



=<JJ+ i V = -[G-^ + Q R {k,u)] 



-i 



(5-11) 
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The imaginary part of Eq. (5-10) is 
lmU R (p,uj) 

= 2 / (0 ( tanh |f " coth ImS ^ + «. « + 1°)^G R (q, q°). (5-12) 

Using Eq. (5-10), the retarded Green function, Eq. (5-2), can be written 

From the rotational and parity invariances, we find that the self-energy has the 
same Dirac matrix structure as the spectral function, 

S R (k,u) = U (k,uj)-f - U v (k,uj)k -7 + E s (k,co), (5-14) 

with S = (l/4)Tr[E R i°], E v = (l/4)Tr[U R k • 7] and E s = (l/A)Ti[E R ]. From 
Eq. (5-10), one can easily check that U s vanishes in the chiral limit, which also 
means that p s = in this case. 

It is useful to decompose £ R (k, lo) and A(k, lo) into positive and negative energy 
parts using the projection operators A T (k) = (1 ± 7% • fe)/2. We obtain 

U R (k,oo) =7°(r_(fc,o;)yl_(fc) + S+(k,u)A+(k)), 
A(k,u) = {p-(k,u)A-(k) + p + (k,Lo)A + (k)) 7 °, (5-15) 



where 



Z T {k,uj) = ^Tr[Z R (k,Lo)A T (k)-f°] = S (k,co) T S v (k,u), (5-16) 
p T (k,u) = ^Tr^^h ^)] = ilm 



2 +v " vr lo + p^k - S T (k,uo) 

= p (k,tu) ± p v (k,Lo) (5-17) 

are the self-energies and spectral functions for positive and negative energies, respec- 
tively. Using Eq. (5-16), the retarded Green function is written as 

G R (k U ) = A - (kh ° + Mkh ° (5-18) 

' u + fj,- k- S-(k,u) +irj lo + p + k - U + (k, lo) + irj 

- + W,. (5.19) 



i?_(fc, to) + it] R + (k,Lo) + ir] 
In the second equality, we have defined 

R T (k,u) = lj + pTk - U T {k,Lo), (5-20) 
for later convenience. The equations 

Re R T (k,co T ) = 0, (5-21) 
give the dispersion relations for quarks and anti-quarks lo = L0 T (k), respectively. 



24 



M. Kitazawa, T. Koide, T. Kunihiro and Y. Nemoto 



There are several possible choices of the four-momentum cutoff for the integral in 
Eq. (5-10). To compute this integral, we first calculate the imaginary part Eq. (5T2); 
substituting the explicit formula for Im G R (k,u>) to Eq. (5T2), the q° integral is 
removed by the delta function in Im G R (k,uj) (see, Eqs. (C-10) and (CT1)). Then, 
Im S R (k,uj) is calculated using the usual three-momentum cutoff A. The real part 
is then evaluated from Im E R (k,uj) through the dispersion relation. To calculate 
Re U R (k,uj) with the dispersion relation, we again need to introduce a cutoff. Here, 
we assume that the cutoff for the dispersion relation is the same as the momentum 
cutoff A, and we then have 



We have checked numerically that the quark spectral function and the DOS near the 
Fermi energy are hardly affected by the choice of the cutoff in Eq. (5-22). 

5.2. The quasiparticle picture of quarks 

Now we present the numerical results and discuss the properties of the quarks 
above T c . To investigate the quasiparticle picture of the quarks, it is useful to consider 
the behavior of the dispersion relation uj = ujzp(k) and the spectral function po(k, uj). 
Here, co^(k) is defined above; see Eq.(5-21). A remark is in order here. Because 
uizf(k) is only the solution to the real part but not the whole part of R^(k,uj^), 
uJzf{k) may not correspond to the peak position of the spectral function and does 
not represent physical excitations when the imaginary part of the Green function is 
large. 

We plot the the dispersion relation of the positive energy quarks uj = u- (k) for 
H = 400 MeV and e = 0.01 in the left panel of Fig. 13. One sees from the figure 
that the dispersion relation deviates from that for free quarks, 10 = k — fi, (dotted 
line) and exhibits a rapid increase around the Fermi energy, uj = 0. We also show 
(duj(k) / d\k\)~ l in the right panel of Fig. 13. This quantity is proportional to the 
DOS, provided that the imaginary part of the self-energy is ignored. As shown in 
the figure, there appears a pronounced minimum of (9cj(fc)/5|fc|) _1 near the Fermi 
momentum, kp = fJ,, reflecting the rapid increase of uj-(k). Although this rapid 
increase may not be physical because there exists a large imaginary part of the 
self-energy in this kinematical region, it may suggest that the gap-like structure is 
induced in the dispersion relation as a precursor to CSC. We will see that this is the 
case from the behavior of the spectral function. 

We show the spectral function po(k,uj) for fi = 400 MeV with e = 0.01 and 
0.2 in Fig. 14. One can see two families of peaks in both figures, near uj = k — p 
and uj = —k — fi, which correspond to the quasiparticle peaks of the quarks and 
anti-quarks, respectively. A notable point is that the quasiparticle peak has a clear 
depression around the Fermi energy, uj = 0, which represents the enhancement of 
the decay rate of the quasiparticles around the Fermi energy. Thus the depression 
becomes more pronounced as e becomes smaller. This behavior is quite different from 
that of the conventional Fermi liquids, for which the lifetime of the quasiparticles 
becomes longer as uj approaches the Fermi energy, uj = 0. 




(5-22) 
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320 360 400 440 480 320 360 400 440 480 

k [MeV] k [MeV] 

Fig. 13. The dispersion relation for the quarks, co = u)-(k) in the left panel, and (<9u>(fc)/<9|fc|)~ 1 
in the right panel with p = 400 MeV and e = (T — T c )/T c = 0.01. Near the Fermi momentum 
\if — H, we see a rapid increase of u_(fc) and a pronounced minimum of (dcu(k)/d\k\)^ 1 around 
k — kF, which means that a gap- like structure is induced in the dispersion relation as a precursor 
to CSC. 



p[MeV 1 ] p[MeV 1 ] 




Fig. 14. The spectral function po at fi = 400 MeV and e = 0.01 and 0.2. The peaks near ui — k — /i 
and lo = — k — fi correspond to the quark and anti-quark quasiparticles, respectively. Note that 
there is a depression around u> = 0, which is responsible for the pseudogap formation. 

To understand these types of non- Fermi liquid behavior of o;_(fc) and pQ(k : uj) 
near to = 0, we investigate the quark self-energy E R (k,io). Before giving the numeri- 
cal results, we note that the spectral function of the negative energy p+(k, uj) always 
takes small values around oj = 0, because the real part of R + (k,u>), Re R + (k,oo) ~ 
uj + /j, + k, does not become smaller around u = 0. Thus, po(k,uj) near the Fermi 
energy is reproduced by p_(fe,cj) alone: 

p-(k,u) + p+(k,u>) 1 
Po(k,uj) = ~-p_(k,uj). (5-23) 

Therefore, it is sufficient to consider the self-energy of the positive energy U^(k,uj) 
in order to understand the non-Fermi liquid behavior of uj^(k) and po(k,uj). 

In Fig. 15, we plot the real and imaginary parts of U-(k,u})\k=k F for p = 400 



26 



M. Kitazawa, T. Koide, T. Kunihiro and Y. Nemoto 



10 




-200 -100 100 200 

co [MeV] 

Fig. 15. The self-energy E_ for k — k F with /i = 400 MeV and e = 0.01. We observe a peak in 
ImS- and a rapid increase of Re£- at the same lu = 0. 




C0h,kh 
• • C0s,ks 



Fig. 16. The decay processes of the quarks described in the non-selfconsistent T-matrix approxi- 
mation. 



MeV and e = 0.01. One finds a peak structure in |Im S-(k, uj)\ and a rapid increase 
in Re U-(k,u>) near u = 0. These phenomena are closely related to the non-Fermi 
liquid behavior of po(k,cu) and uj-(k). In fact, the former means that the decay rate 
of the quasiparticles is enhanced around the Fermi energy, and the latter implies an 
increase of uj-(k) around uj = 0. 

To understand the origin of the characteristic behavior of I7_ (k, u>) around oj = 0, 
let us consider the decay mechanism of a quark. In our calculation, a quark decays 
only through the process depicted in Fig. 16, where an incident quark takes up 
another quark to make a hole "h" and a diquark "(qq)" state; 

q^h+(qq). (5-24) 

This process is enhanced when the diquark pair is in a collective state, provided that 
the energy-momentum matching is satisfied, i.e. 

(u,k) = (uj h ,k h ) + (uj s ,k s ), (5-25) 

where (u, k), (ujh,kh) and (iv s ,k s ) denote the energy-momentum of the decaying 
particle, the hole, and the collective mode composed of the diquark pair, respectively. 
We have seen in § 3 that the diquark pair field makes a collective soft mode and the 
pair fluctuations are strong near u = k = as T approaches T c . Thus one sees that 
the decay process Eq. (5-24) is enhanced for (cj s ,fc s ) ~ (0,0) when T is close to T c . 
On the other hand, the hole energy should become u>h = M — \kh\, since the hole 
in Eq. (5-24) is on-shell. Combining these two conditions, we find that the decay 
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Fig. 17. The imaginary and real parts of the quark self-energy E-(k,ui) with /i = 400 MeV and 
e = 0.01. One observes a peak in llml?- 1 and a rapid increase of ReS- near to — fi — k. 

process Eq. (5-24) is most enhanced when 

(w,fc) ~ (w h ,k h ) = (ji- \k\,k), (5-26) 

where we have used the momentum conservation in the last equality. We show 
that this is indeed the case in the left panel of Fig. 17, where we show the energy- 
momentum dependence of Im U_ (k, uj) for fi = 400 MeV and e = 0.01. There we can 
clearly see peaks along the line uj = fi — \k\, as expected. At the Fermi momentum 
k = kF, this peak corresponds to uj = 0, as shown in Fig. 15. We also show in the 
right panel of Fig. 17 the energy-momentum dependence of Re E-{k, uj) for the same 
\i and e. There a rapid increase in Re E-(k,u) is seen along the line uj = fi— \k\, at 
which |Im E-(k,u)\ has a peak structure, in accordance with the behavior expected 
from the dispersion relation Eq. (5-22). 

Substituting po into Eq. (5-4), we obtain the DOS, N(u), of the quarks including 
the effect of the pair fluctuations. In Fig. 18, we show the DOS at fx = 400 MeV for 
e = 0.01, 0.02, 0.05 and 0.2. The DOS of the free quarks 

, . 2NrN c , N . 
N (uj) = ^^(uj-fi) 2 , (5-27) 

is also shown by the thin dotted curves for comparison. We find that there appears 
a clear depression in the DOS around the Fermi energy for e = 0.01. 27 ) This is the 
pseudogap of CSC. This pseudogap survives up to e ~ 0.05. Therefore, the pseudo- 
gap certainly does appear as an effect of the pair fluctuations. The appearance of the 
pseudogap in the quark DOS can be naturally understood as a manifestation of the 
non-Fermi liquid behavior in po(k,uj) and uj = uj-(k), because both the pronounced 
minimum of the quasiparticle peak and the rapid increase of the dispersion relation 
around uj = indicate a decrease of the DOS near uj = 0. Although we did not 
calculate the DOS for smaller e, say, e < 0.01, because of numerical difficulties, it 
seems that the gap-like structure grows as e is lowered. In such a region, however, 
the non-linear effects of the pair fluctuations become significant, and a systematic in- 
corporation of theses nonlinear effects using, for instance, the renormalization group 
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Fig. 18. The density of states for n = 400 MeV and various e. The thin dotted curve represents 
that of free quarks. A clear pseudogap structure is seen, which survives up to e ~ 0.05. 
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Fig. 19. The chemical potential dependence of the DOS. Here, the quantities plotted is the ratio 
of the DOS for the system under study to that for free quarks, No{u>). The pseudogap becomes 
more pronounced as /i increases. 



analysis, should be made to obtain the correct behavior of the DOS. Such a treatment 
is, however, beyond the scope of the present work. 

To elucidate the chemical potential dependence of the DOS, we show the DOS 
for fi = 350, 400, 500 MeV in Fig. 19, where the values actually plotted are the ratios 
N(uj)/No(uj) of the DOS for the system under study to that of free quarks given 
in Eq. (5-27). From Fig. 19, we can conclude that there is a "pseudogap region" 
within the QGP phase above T c up to T* = (1.05 - 1.1)T C at intermediate densities. 
We also find that the pseudogap becomes more significant as fx increases. This is 
because the Cooper instability becomes stronger as /j, increases owing to the larger 
Fermi surface. 

A comment is in order. In Ref. 65), the pseudogap in low density nuclear matter 
is investigated, and it is found that a clear pseudogap is seen up to e « 0.0025, 
which is more than one order of magnitude smaller (in units of e) than our result. 
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This is simply a reflection of the strong coupling nature of QCD in the intermediate 
density region. Our result obtained for a three-dimensional system reveals that a 
considerable pseudogap can be formed without a low-dimensionality effect, as in 
the HTSC, and that the pseudogap phenomena may be universal in strong coupling 
sup er conductivity . 64 ) 

§6. Summary and concluding remarks 

In the present paper, we have examined precursory fluctuations of the diquark- 
pair field above the critical temperature T c for two-flavor color superconductivity 
(CSC) in heated quark matter at moderate density A Nambu-Jona-Lasinio-type 
model was adopted as an effective theory embodying the strong coupling nature of 
QCD, which becomes more relevant for quark matter at such a relatively low den- 
sity. The detailed formulation given here is based on linear-response theory, and a 
detailed account of the calculational procedure has also been given in the imaginary- 
time formalism. We showed that the pair fluctuations develop a collective mode 
whose complex frequency approaches the origin of the complex energy plane as T 
decreases toward T c ; i.e., the pairing fluctuations form the soft mode of CSC. More- 
over, it was shown that the pairing fluctuations are quite strong even well above T c , 
owing to the strong-coupling nature of the dynamics, in comparison with usual metal 
superconductors. We have presented an extensive investigation of the properties of 
the precursory pair field with finite momenta. We calculated the spectral function 
of the pair fluctuations p(k,u>) as a function of the momentum k and energy to to 
clarify the spatial and temporal behavior of the pair fluctuations. We also presented 
the behavior of the dynamical structure factor. It was found that the collective mode 
with a shorter wavelength has a larger decay width and smaller strength. 

We examined the effects of the precursory pair fluctuations on the specific heat 
and the quark spectrum at values of T close to but above T c . It was shown that 
c v increases rapidly when T is lowered toward T c , and it diverges at T = T c , along 
with singular growth of the pair fluctuations, in accordance with the general theory 
of second-order phase transitions. 

The single-quark Green function was also calculated for the first time by incor- 
porating the effects of the diquark-pair fluctuations in the T-matrix approximation. 
We showed that the quarks behave like a typical non-Fermi liquid, owing to the 
soft mode. In particular, the quarks have a larger width near the Fermi surface as 
a result of the interaction with the pairing soft mode. This leads to a pseudogap 
in the density of states (DOS) of the quarks in the vicinity of the critical point. 
The chemical potential dependence of the quark spectrum and the DOS were also 
studied. Our results suggest that the heated quark matter in the vicinity of CSC 
phase transition at moderate densities is not a simple Fermi liquid, due to the strong 
fluctuations of the pair field. This may invalidate the mean-field approximation. It 
is also to be noted that a pseudogap can form even in (relativistic) three-dimensional 
systems without the need for low-dimensionality, as has been suggested for HTSC. 
This leads us to conjecture that the formation of a pseudogap may be a universal 
phenomenon for strongly correlated matter irrespective of the spatial dimensionality 
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A 




Fig. 20. The diagrams that contribute to the photon (or gluon) self-energy, representing the 
Aslamazov-Larkin term (left) and the Maki-Thompson term (right). 

of the system. 

What observables can be most strongly affected by the precursory fluctuations 
of the diquark-pair field? It is known in the context of the condensed matter physics 
that pair fluctuations above T c affect several transport coefficients, including the 
electric conductivity (EC) and the phonon absorption coefficient, as well as the spe- 
cific heat. 66 ) A large excess of electric conductivity seen near but above T c is known 
as paraconductivity. 66 ) They are attributed to pairing fluctuations in the normal 
phase: Typical microscopic mechanisms that can give rise to such an anomalously 
large conductivity are those resulting from the so-called Aslamazov-Larkin and Maki- 
Thompson terms. 67 -* Those are depicted in Fig. 20. The dotted lines in the figure 
denote the gauge field, i.e., the electro- magnetic field (or the photon) in this case. A 
direct analogy of electric superconductivity to CSC may be realized by replacing the 
photon field by the gluon field in Fig. 20, leading to color-paraconductivity above T c in 
our case. However, unfortunately it would be difficult detect the color-conductivity 
directly in experiment or observation. Nevertheless, note, however, that the external 
photon field can couple to the diquark-pair fluctuations of CSC. This coupling is 
precisely depicted by the same diagram, Fig. 20, with the pair field interpreted as 
being composed of colored quarks; in other words, the photon self-energy in quark 
matter at T > T c is strongly modified due to the fluctuating diquark pair field. This 
is interesting, because modifications of the photon self-energy in heated quark matter 
near T c should bring about an enhancement of dilepton emission from the system, 
which may carry some information concerning the fluctuations of the diquark-pair 
field in heated quark matter, created, say, in the intermediate stage of heavy-ion 
collisions. It is also conceivable that the electric conductivity in quark matter will 
show anomalous enhancement at T near T c in the normal phase due to fluctuations 
of the colored pair-field through the same process. 

The precursory phenomena may also affect the cooling process of compact stars. 
The temperature of newborn compact stars just after a supernovae can exceed 40 
MeV. Therefore, dense matter in the interior of compact stars may undergo CSC 
phase transition. However, it may be the case that the /3-equilibrium condition is 
not completely satisfied in proto-compact stars just after supernova explosion, in 
contrast to the situation in the interiors of cold compact stars, as mentioned in the 
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Introduction. Furthermore, a soft mode necessarily appears if the phase transition 
is second order, irrespective of the pattern of CSC. Therefore, the results presented 
in the preceding sections can apply to the new-born compact stars. If it is the case, 
then the precursory phenomena could affect their cooling process. Indeed, we have 
shown that the specific heat is enhanced through the precursory diquark soft mode 
as T approaches T c from above. Moreover, the neutrino mean-free path should be 
affected and become shortened through the scattering with the soft mode near T c . 

In the present work, we have employed the non-self-consistent T-matrix approx- 
imation. As another, more complicated approximation, there is the self-consistent 
T-matrix approximation, in which all free fermion propagators are replaced by re- 
summed fermion propagators. However, the self-consistent approach does not nec- 
essarily yield better approximation. As shown in Ref. 68), the vertex corrections to 
the self-energy ignored in the self-consistent approximation cancel each other, and 
only the lowest-order term survives. As shown in Ref. 68) again, if one takes into 
account the vertex corrections as well as the self-consistency to the self-energy, their 
contributions approximately cancel each other, and correspondingly the results are 
rather close to those obtained with the non-self-consistent approximation. 

In this paper, we have used a fixed diquark coupling Gc to examine the T and 
fi dependences of the quark spectrum and the DOS. An investigation of the quark 
spectrum with a larger Gc may provide a deeper understanding of the underlying 
physics for the formation of the pseudogap, for instance, in terms of the so called 
resonance scattering. 30 )' 64 )' 69 ) 
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In this appendix, we calculate the quark-quark polarization function Q{k,u n ) 
given in Eq. (3-15), which reads 



Appendix A 

Calculation of Q{k,u) 




Tr [i75T2A 2 CC/o(fc ~P,v n - w^i^T-i^CQ^ (p, cj. 



rn 



)] 




rn 



3 Trc [g (k -p,u n - u m )Go(p, Lo m )] , (A-l) 
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where TYp denotes the trace over the Dirac spinor. Here, Qo{k,u n ) is the free 
Matsubara Green function 

g (k,u n ) ee 1 = £ , f s7 [ , , , (A-2) 

where A T = (1 ± 7°7 • k)/2 are the projection operators onto positive and negative 
energies, respectively. Substituting Eq. (A-2) into Eq. (A-l), we obtain 

J '''' — r \k — p\j° + s(k — p) • 7 |p|7° + tp ■ 7 



m y 1 s,t=± 

1 



2|fc-p| 2|p| 



iv n - iuj m + n + s\k — p\ iiv m + fj, + i|p| 
= 2iVfA^-nV / rf3 P l fc ~PlN -st(fc-p) -p/ + (-s|fc-p|) -/~(t|p|) 
" s^±^ (27r)3 iv n + 2n + s\k-p\+t\p\ 

N ( N _n v / d3p f (i fc -pi + ipi) 2 - fc2 -p|)-rHp|) 

J s tl^( 2vr ) 3 l M/„ + 2 M + a(|fc-p| + |p|) 

(|fe - P\ ~ \P\? ~ k 2 f + (s\k - p|) - /-(- S |p 



— p||p| iu n + 2fi + s(\k - p\ - \p\) 
N f (N c -l)^ f^_fr, , , 2 , 2l / + (- S ei)-/-( S e 2 ) 



(A-3) 



/ de 1 de 2 \ [(ei + e 2 ) 2 - fc 2 l 



- [(ei - e 2 ) 2 - fc 2 ] /+( - Sei) - r( - Se2 H . (A4) 
LV ; J + 2/x - s(ei -e 2 ) J V ; 

Using Eq. (A-3), Eq. (3-18) can be obtained with a simple manipulation. In the last 
equality, we have introduced the new variables e\ = \k — p\ and e 2 = |p|, which are 
such that 

I T^zTi ~7TT = /o \ti I de i de 2, (A-5) 

J (2vr)^ |k-p||p| (2nyk J s 

with S being the domain of integration shown in Fig. 21 a). 

Performing the analytic continuation Q R (k,u>) = Q(k,v n )\i Un=U) +i v and using 
the fact that the domain S is invariant under the exchange of e± and e 2 , we have 

X ([(.,+ erf" -flf ) 

[ \lo + 2ix - e\ - e 2 + if] uj + 2ix + ei + e 2 + vq J 

+2[( ei - e 2 ) 2 - fc 2 ] /'( e i)-/ + ( 62 ) 1 (A . 6) 

N f (N c -l) f dE 1 dE 2 
~ (2vr) 2 A; 
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Fig. 21. The domains of integration in Eqs. (A-4) and (A-7). 



x (E\ - k 2 ) 



tanh ^i+g-^ 



to + 2/x — + ir] uj + 2fi + Ei + ir] 



f-( Ei+E 2 \ _ r+( E 1 -E 2 \ 

+2(E 2 - k 2 ) 1 { \ ' I { 2 ' 



(A-7) 



In the second equality, the new variables E\ = e± + e 2 and E 2 = e± — e 2 have been 
introduced for later convenience to obtain the imaginary part of Q R (k,oo). 

As explained in the main text(§3), we first calculate the imaginary part Q with- 
out introducing a cutoff, and then we define the real part using the dispersion relation 
with a momentum cutoff introduced. From Eq. (A-7), we have for the imaginary part 
of Q R , 



ImQ R (k,uj) 

N f (N c -l) 



= — 7T- 



87T 2 k 



oo />k 

dE-i I dE 2 
k J-k 



x {{E\- k 2 ) (tanh Ei+ ^- 2 » 5(lu + 2fi- E ± ) - tanh El+ g +2/J 5{u + 2fi + £i)) 
+2{E 2 - k 2 ){f(^ - M ) - /(^ + + 2^ - £ 2 )} 

= - ^^"^ {^(k + M ~ k )j dE 2 [{to + 2/x) 2 - /c 2 ] tanh ^ 

+20(-|w + 2/x| + k) dE 1 [(u + 2^) 2 - k 2 ] {/(^) - /(^)}| 

Nf(N c — 1)T r n n2 i21 r,. n cosh(u; + /c)/4T 

= — — ^ + 2 A*) - k 2 ] < 9(\u + 2d - fc) log — -j -^j— 

2nk LV ^' J \ Vl ^ ; cosh (w - fc)/4T 



+0(-|w + 2^| + fc)log 



l + e -(^+fc)/2T | 

1 + e -(-u+k)/2T j 



N f (N c -l)T 



[(oj + 2fi) 2 - k 2 } 



X <^ loi 



2irk 

cosh(w + fe)/4T w 



cosh(w - fc)/4T 2T 



- — 0(-\u + 2n\ +k) } . 



(A-8) 
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As stated above, we next define the real part of Q R through the dispersion 
relation with Eq. (A-8) for the imaginary part: 

ReQ R (k,uj) = -h [duj MQ ^f . (A-9) 

TV J UJ — CO 

We introduce a cutoff so that the integration of Eq. (A-9) becomes consistent with 
the Thouless criterion Eq. (3-21). With vanishing momentum, Eq. (A-8) becomes 

ImQ R (k -> 0, to) = - Nf ^ N 4 c - ^ (to + 2fi) 2 tanh + 0{k 2 ). (A-10) 

47r 4T 

Substituting this equation into Eq. (A-9), we have 

n „, . Nt(N c -l) f (uj + 2/j,) 2 tanh £L 

RcQ R (0,0) = n A \ / duj- (A-11) 

4n z J uj 

In comparison with Eq. (2-9), we find that the cutoff should be introduced as in 
Eq. (3-23). 

Appendix B 

Analytic Properties of Q(k,uj) 



In order to find the pole of the response function D R (k,u>), we have to derive 
the analytic form of Q R (k,uj) in the lower-half complex energy plane, as given in 
Eq. (3-26). 

We here restate the definitions of Q(k, z), Q R (k, z) and Q A (k, z) in the complex- 
energy plane C 

Q(k,z) = Q(k,is n )\ lUn ^ z , (B-l) 
Q R (k,z) = Q R (k,uj)\^ z , (B-2) 
Q A (k,z) = Q A (k,uj)\^ z . (B-3) 

Because Q(k,z) has a cut along the real axis, Q R (k,z) (Q A (k,z)) coincides with 
Q(k, z) only in the upper (lower) half plane C + (C~): 

Q{k,z) = Q R (k,z) for z G C + , 

Q(k,z) = Q A {k,z) for zeC~. (B-4) 

In this appendix, we wish to find the analytic form of Q R (k,z) in C~. To 
accomplish this, we first evaluate Q R (k,z) just below the real axis, at z = uj — ir), 
and then determine Q R (k,z) in C~ by analytic continuation. Here, we note that 
the discontinuity of Q(k,z) at the cut along the real axis is characterized by the 
discontinuity of the imaginary part, that is, 

ReQ R (k, uj + irj) = ReQ(k, uj + irf) = ReQ(k, uj - irf) = ReQ A (k, to -irf), (B-5) 
lm.Q R {k, uo + irf) = ImQ(k, to + irf) = —lmQ(k, uj — irf) = — lmQ A (k, uj — ir]) 

= I(k,u), (B-6) 
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(a) 


1 k 


(b) 


> 





V 



Fig. 22. The analytic continuation of Q R (k,uj) for a vanishing momentum. 

with lo being a real variable. We find from Eq. (3-22) that cuts of I(k, z) = 
I(k 

)^)|cj >z ©xist git Imz — 2(2ti -|- \^ttT. HencG, this function is analytic in Cj = 

{z G C| - 2ivrT < Imz < 2ivrT}. Using Eqs.(B-5) and (B-6), we find 

Q R (k, lo) - il(k, lo) = Q R {k, Lo + irj) - ilmQ R (k, lo + in) 
= ReQ R (k, uj + if]) = ReQ A (k, lo - irj) 
= Q A (k, lo — irj) — ilmQ A (k, to — ir]), 
= Q(k,oj -irj) +iI(k,uo), (B-7) 



which leads to 



Q R (k, oo - irj) = Q R (k, to + irj) 

= Q(k, lo — irj) + 2il(k, lo). 



(B-8) 



In the first equality here, we have used the fact that Q R (k, z) is a continuous function 
on the real axis, and Eq. (B-7) is applied to derive the second equality. Equation 
(B-8) implies that we can express Q R (k,z) just below the real axis in terms of the 
functions Q(k,z) and I(k,z), whose analytic forms in Imz < are already known. 
Because both Q(k,z) and I(k,z) are analytic functions in C~ PI Cj, the analytic 
form of Q R in this domain can be obtained straightforwardly from the uniqueness 
theorem of the analytic continuation as 



Q R (k,z) =Q(k,z) + 2iI(k,z) 



for z G C~ n Cj. 



(B-9) 



Although Eq. (B-9) is valid only for Imz > —2ttT, this region is sufficient to find the 
pole. 

For vanishing momentum, analytic continuation into the lower half plane is 
understood as changing the contour of the momentum integration as shown in Fig. 22. 
In this case, the term originating from the pole integration is added in the lower half 
plane, 8 ) and we obtain 



Q R {0, z G C~) 



2N f (N c -l) 



71 



|k|-M 



r P [ A k 2 dk- tanh ~ 

J-A Z 



2(\k\-ji) 



- 2 2m ^ + 2 



Z 2 Z 

1 tanh 



4T J 



(B-10) 



Equation (B-10) coincides, of course, with the k —> limit of Eq. (B-9). 
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Appendix C 

Calculation of U 

To carry out the summation of the Matsubara frequency uo' m in Eq. (5-7), it is 
sufficient to consider the case uj n > 0. Because both S(k,z) and Go(k,z) have a 
cut along the real axis, the integrand of Eq. (5-7) has two cuts along the horizontal 
lines Imz = and Imz = —iu> n in the complex energy plane as shown in Fig. 23. 
Therefore, we must to separate the Matsubara frequency into four parts to carry out 
the Matsubara sum in Eq. (5-7): 

S(p,uj n ) = 8T^2 I ^-j3^(p + q,Wn + w m )e (9,w m ) 

R (p + q,ioj n + q )G R (q,q°) 

R (p + q,iLO n + q°)G A (q,q ) 

(p + q,iuj n +q )G A (q,q ) 

A (p + qi iuj n + q )G£(q,q ), (C-1) 

where the contours C\ - C4 are shown in the left panel of Fig. 23. The definitions of 
the retarded and advanced T-matrices E R and S A are 

S(p, z) = S R (p, z) for z G C + , 

~(p, z) = S A (p, z) for z G C-, (C-2) 

Next, we change the contours C\ - C4, as shown in Fig. 23. Then the contribution 
from the integral along the large circle vanishes, and only the integrals along the 



J Cl 2m 2T J 

+ i f ^tan„4/ 
J Cl 2m 2T J 

+ i f ^ t anh£/ 
J Ca 2m 2T J 

J Ct 2m 2T J 
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horizontal lines Imoj = and Imcj = — iu n remain. Thus, we have 
~, . f d 3 q f°° dq° , q° 

xS R (p + q, iu n + q°) {G R (q, q°) - G A (q, q )} 
f d 3 q [<*>-*»» dq° , q° 

X {S R (p + q, ioj m + q°) - ~ A (p + q, iuj m + q )} G A (q, q°) 
= 2 j ^{^^ R (P + <l^n + q )lraG R (q, q ) 

+ tanh l^pLl m ~ R (p + q, q°)G A (q, q° - iu n ) } , (C-3) 

where / d 4 g/(2vr) 4 = / dq°/(2ir) J d 3 g/(2vr) 3 , and the relations 

Im~ R = (~ R - ~ A ) /2i, lmG R = (G R - G A )/2i (C-4) 

have been used in Eq. (C-3). It is known that the singularity corresponding to 
the integral along C3 does not give rise to any problem. 55 ) Then, bearing in mind 
that tanh - l££ n = coth 7^, and employing the analytic continuation E R {k,u) = 
, ^n)\iu!„=u>+ir]-, we obtain the self-energy in real time, 

E R (k,u) = 2 J ^|tanh^^(p + q^ + ^)ImG ^,g°) 

+ coth ^lmS R (p + q, q°)G A (q, q° - u) \ . (C-5) 

To simplify Eq. (C-5) further, we substitute the explicit form of the free Green 
function, 



G { RA} ( q ) = fi A - {q) , + „ — , , » 
q° + H±ir]-\q\ q° + fi ± ir] + \q\ 

and its imaginary part, 

lmG { R ' A} (q) = Tvr {A-{q)S(q° + fi - \q\) + A + (q)5(q° + » + \q\)) , (C-7) 
into Eq. (C-5). Then, Eq. (C-5) becomes 

U (p,u J ) = ^TT{ 1 U(p,oj)} 

= A [A, I P lm f(P + ^ [tanh JgLz£- coth ^] 
27 (2vr) 3 7 27ru;-p -|g|+/x- V 2T 2T J 

^ lrnS R (p + q,u) -\q\ - M _ a, 

27 (2vr) 3 7 2tt a; + |q| + ^ - iT7 L 2T 2T J ' 

(C-8) 
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4 

f d 3 q f du Im S R (p + q, uj) A Igl — a . u , 

= - / 7—77 / n — i i (g-p)tanh^ — - coth — 

J (27r) 3 7 2iruj-p°- \q\+(i-iri yH y>V 2T 2T J 

+ [d*q [ p hnSy + q, U ) -\q\ - „ _ a, 

7 (2tt) 3 7 2vru;-p + jql + ^-ir?^ 2T 2T J 

(C-9) 

The imaginary parts of these quantities are found to be 
lmU (p,p°) 

= -\f |^Im =*(p + + |«| - /,)[tanh ^ - coth H^fa] 

"I / ^ Im ^(P + 9,P°-|Q|-^)[tanh^-eoth^fa], (C-10) 
Im U v (p,p°) 

= ~\j T^Im + q,p° + |q| - »)(q ■ p)[tanh ^ _ coth Alto] 

+ I / (C^ Im ^ + " N " ^ " ^ )[tanh ^ " C ° th £!Z ^ ] - 

(C-ll) 
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